switch dataswitch
    case 'T2eff7T'
        temp = 1./niftiread('R2star-7T.nii.gz');
    case 'R2eff7T'
        temp = niftiread('R2star-7T.nii.gz');
    case 'F'
        temp = niftiread('qMT-BSP-Parameter-F.nii.gz');
    case 'R1obs7T'
        temp = niftiread('R1-7T.nii.gz');
    case 'T1obs7T'
        temp = 1./niftiread('R1-7T.nii.gz');
    case 'R1obs3T'
        temp = niftiread('R1-3T.nii.gz');
    case 'T1obs3T'
        temp = 1./niftiread('R1-3T.nii.gz');
    case 'Amy'
        temp = niftiread('Histology-Gallyas.nii.gz');
        tempmask = zeros(size(temp));
        tempmask(:,:,50) = mask;
        temp = temp./max(temp(tempmask==1));
    case 'AFe'
        temp = niftiread('Histology-Perls.nii.gz');
        tempmask = zeros(size(temp));
        tempmask(:,:,50) = mask;
        temp = temp./max(temp(tempmask==1));
    case 'RM0a'
        temp = niftiread('qMT-BSP-Parameter-RM0a.nii.gz');
    case 'R1a'
        temp = niftiread('qMT-BSP-Parameter-R1a.nii.gz');
    case 'R2a'
        temp = niftiread('qMT-BSP-Parameter-R2a.nii.gz')./10^3;
    case 'R2b'
        temp = niftiread('qMT-BSP-Parameter-R2b.nii.gz')./10^6;
    case 'T1a'
        temp = 1./niftiread('qMT-BSP-Parameter-R1a.nii.gz');
    case 'T2b'
        temp = 1./niftiread('qMT-BSP-Parameter-R2b.nii.gz').*10^6;
    case 'T2a'
        temp = 1./niftiread('qMT-BSP-Parameter-R2a.nii.gz').*10^3;
end